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We study some exact and approximate solutions of Einstein's equations that can be used to 
describe the gravitational field of astrophysical compact objects in the limiting case of slow rotation 
and slight deformation. First, we show that none of the standard models obtained by using Fock's 
method can be used as an interior source for the approximate exterior Kerr solution. We then use 
Fock's method to derive a generalized interior solution, and also an exterior solution that turns out 
to be equivalent to the exterior Hartle-Thorne approximate solution that, in turn, is equivalent to 
an approximate limiting case of the exact Quevedo-Mashhoon solution. As a result we obtain an 

£Nj ■ analytic approximate solution that describes the interior and exterior gravitational field of a slowly 

t-H 1 rotating and slightly deformed astrophysical object. 
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, ■ Keywords: Hartle-Thorne metrics, Fock metric, Kerr metric, Quevedo-Mashhoon metric 

hi, 

I. INTRODUCTION 

qh In astrophysics, the term compact object is used to refer to objects which are small for their mass. In a wider sense, 
I ■ the class of compact objects is often defined to contain collectively planet-like objects, white dwarfs, neutron stars, 
£^ other exotic dense stars, and black holes. It is well known that Newtonian theory of gravitation provides an adequate 
i— — i. description of the gravitational field of conventional astrophysical objects. However, the discovery of exotic compact 
objects such as quasars and pulsars together with the possibility of continued gravitational collapse to a black hole 
' . points to the importance of relativistic gravitation in astrophysics. Moreover, advances in space exploration and the 
■ development of modern measuring techniques have made it necessary to take relativistic effects into account even 
\ in the Solar system. Probably the simplest way to study the relativistic gravitational field of astrophysical compact 
f^S ■ objects is by expressing it in terms of their multipole moments, in close analogy with the Newtonian theory, taking 
C*~) ' into account the rotation and the internal structure of the source. 

^ | . In this context, the first exterior solution with only a monopole moment was discovered by Schwarzschild 

1 soon after the formulation of Einstein's theory of gravity. In 1917, Weyl @ showed that the problem of finding 
CNJ ■ static axisymmctric vacuum solutions can generically be reduced to a single linear differential equation whose general 
t-H ■ solution can be represented as an infinite series. The explicit form of this solution resembles the corresponding solution 
\ in Newtonian's gravity, indicating the possibility of describing the gravitational field by means of multipole moments. 
In 1918, Lense and Thirring Q discovered an approximate exterior solution which, apart from the mass monopole, 
contains an additional parameter that can be interpreted as representing the angular momentum of the massive body. 
From this solution it became clear that, in Einstein's relativistic theory, rotation generates a gravitational field that 
leads to the dragging of inertial frames (Lense-Thirring effect). This is the so-called gravitomagnctic field which is 
of especial importance in the case of rapidly rotating compact objects. The case of a static axisymmctric solution 
with monopole and quadrupole moment was analyzed in 1959 by Erez and Rosen [H by using spheroidal coordinates 
which arc specially adapted to describe the gravitational field of non-sphcrically symmetric bodies. The exact exterior 
solution which considers arbitrary values for the angular momentum was found by Kerr [j| only in 1963. The problem 
of finding exact solutions changed dramatically after Ernst [f| discovered in 1968 a new representation of the field 
equations for stationary axisymmetric vacuum solutions. In fact, this new representation was the starting point to 
investigate the Lie symmetries of the field equations. Today, it is known that for this special case the field equations 
are completely intcgrable and solutions can be obtained by using the modern solution generating techniques Q. There 
are several solutions with higher multipole moments j8l4l2l| with very interestingphysical properties. In this work, we 
will analyze a particular class of solutions, derived by Quevedo and Mashhoon [131 ] in 1991, which in the most general 
case contains infinite sets of gravitational and electromagnetic multipole moments. Hereafter this solution will be 
denoted as the QM solution. 

As for the interior gravitational field of compact objects, the situation is more complicated. There exists in the 
literature a reasonable number of interior spherically symmetric solutions fl4| that can be matched with the exterior 
Schwarzschild metric. Nevertheless, a major problem of classical general relativity consists in finding a physically 
reasonable interior solution for the exterior Kerr metric. Although it is possible to match numerically the Kerr 
solution with the interior field of an infinitely tiny rotating disk of dust fl5| , such a solution cannot be used in general 
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to describe astro-physical compact objects. It is now widely believed that the Kerr solution is not appropriate to 
describe the exterior field of rapidly rotating compact objects. Indeed, the Kerr metric takes into account the total 
mass and the angular momentum of the body. However, the quadrupolc moment is an additional characteristic of any 
realistic body which should be considered in order to correctly describe the gravitational field. As a consequence, the 
multipole moments of the field created by a rapidly rotating compact object are different from the multipole moments 
of the Kerr metric [l6j. For this reason a solution with arbitrary sets of multipole moments, such as the QM solution, 
can be used to describe the exterior field of arbitrarily rotating mass distributions. 

In the case of slowly rotating compact objects it is possible to find approximate interior solutions with physically 
meaningful energy-momentum tensors and state equations. Because of its physical importance, in this work we will 
review the Hartle-Thorne [l7l [l8j interior solution which are coupled to an approximate exterior metric. Hereafter 
this solution will be denoted as the HT solution. One of the most important characteristics of this family of solutions 
is that the corresponding equation of state has been constructed using realistic models for the internal structure of 
relativistic stars. Semi-analytical and numerical generalizations of the HT metrics with more sophisticated equations 
of state have been proposed by different authors [l9j]. A comprehensive review of these solutions is given in [lf|. In 
all these cases, however, it is assumed that the multipole moments (quadrupole and octupole) are relatively small and 
that the rotation is slow. 

To study the physical properties of solutions of Einstein's equations, Fock [2(| proposed an alternative method in 
which the parameters entering the exterior metric are derived by using physical models for the internal structure of 
the body. In this manner, the significance of the exterior parameters become more plausible and the possibility appear 
of determining certain aspects of the interior structure of the object by using observations performed in the exterior 
region of the body. Fock's metric in its first-order approximation was recently generalized in 1985 by Abdildin [2"D.|22T| 
to include the case of rotating objects. 

In this work, we review the main exact and approximate metrics which can be used to study the interior and exterior 
gravitational field of compact objects and find the relationships between them. We will show that the exterior HT 
approximate solution is equivalent to a special case of the QM solution in the limit of a slowly rotating slightly 
deformed compact object (first order in the quadrupole and second order in the angular momentum). Moreover, we 
will show that a particular case of the extended Fock metric is equivalent to the approximate exterior HT solution. 
Furthermore, since the exterior HT solution can be matched to the interior HT solution, we conclude that the exterior 
QM solution describes the gravitational field of astrophysical compact bodies and it should be possible to match it 
with an exact interior still unknown solution. 

This paper is organized as follows. In Section [XT] we review the HT solutions and briefly comment on their most 
important properties. In Section Mil we present Fock's extended metric, as first derived by Abdildin [U [22[ in 
harmonic coordinates, and introduce a set of new coordinates which makes it suitable for comparison with other 
exterior metrics. Moreover, we find explicitly the coordinate transformation that relates Fock's extended metric with 
the exterior HT solution. 

In Section HVBI we present a particular case of the QM metric which contains, in addition to the mass and angular 
momentum parameters, an additional parameter related to the mass quadrupole of the source. Here we show explicitly 
that a limiting case of the QM metric contains the HT metric. Finally, Section IVl contains discussions of our results 
and suggestions for further research. 

II. THE HARTLE-THORNE METRICS 

To second order in the angular velocity, the structure of compact objects can be approximately described by the 
mass, angular momentum and quadrupole moment. An important consequence of this approximation is that the 
equilibrium equations reduce to a set of ordinary differential equations. Hartlc and Thorne (l7l . ITsj explored the 
gravitational field of rotating stars in this slow rotation approximation. This formalism can be applied to most 
compact objects including pulsars with millisecond rotational periods, but it shows "large" discrepancies in the case 
of rapidly rotating relativistic objects near the mass-shedding limit [16| . i. e., when the angular velocity of the object 
reaches the angular velocity of a particle in a circular Keplerian orbit at the equator. An additional property of 
this formalism is that it can be used to match an interior solution with an approximate exterior solution. In the 
following subsections we will present the interior and the exterior metrics and introduce notations which will be used 
throughout the paper. 

A. The interior solution 

If a compact object is rotating slowly, the calculation of its equilibrium properties reduces drastically because it 
can be considered as a linear perturbation of an already-known non-rotating configuration. This is the main idea of 
Hartlc's formalism JT^]. To simplify the computation the following conditions are assumed to be satisfied. 

1) Equation of state: The matter in equilibrium configuration is assumed to satisfy a one-parameter equation of state, 
V = V(£), where V is the pressure and £ is the mass-energy density. 

2) Axial and reflection symmetry: The configuration is symmetric with respect to an arbitrary axis which can be 
taken as the rotation axis. Furthermore, the rotating object should be invariant with respect to reflections about a 
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plane perpendicular to the axis of rotation. 

3) Uniform rotation: Only uniformly rotating configurations are considered since it is known that configurations that 
minimize the total mass-energy (e.g., all stable configurations) must rotate uniformly (23|. 

4) Slow rotation: This means that angular velocities H are small enough so that the fractional changes in pressure, 
energy density and gravitational field due to the rotation are all less than unity, i.e. 

where M' is the mass and R' is the radius of the non-rotating configuration. The above condition is equivalent to the 
physical requirement f2 <C c/R'. 

When the equilibrium configuration described above is set into slow rotation, the geometry of space-time around it 
and its interior distribution of stress-energy are changed. With an appropriate choice of coordinates, the perturbed 
geometry is described by 

*» - e* [1 + * + WW, ^ - ^^'t^ln - 2 "' )] ^- (2) 
-R 2 [1 + 2(v 2 - h 2 )P 2 ] [dO 2 + sin 2 9(# - ujdt) 2 ] + 0(D 3 ) 

Here M' is the mass of the non- rotating star, P 2 = P 2 (cosQ) is the Legendre polynomial of second order, u> is the 
angular velocity of the local inertial frame, which is a function of R and is proportional to the star's angular velocity 
f2, and, finally, ho,h 2 , mo, m 2 , v 2 are all functions of R that are proportional to fi 2 . 

In the above coordinate system the fluid inside the star moves with a 4- velocity corresponding to uniform and rigid 
rotation (24| . The contravariant components are 

u* - {gu + Zng^ + rfg^)- 1 ' 2 , = flu' , u R = u e = 0. (3) 

The quantity f2 that appears in the expression for it* is so defined that Q = Q — u) is the angular velocity of the fluid 
relative to the local inertial frame. 

The energy density and the pressure of the fluid are affected by the rotation because it deforms the compact object. 
In the interior of the object at a given (i?, O), in a reference frame that is momentarily moving with the fluid, the 
pressure is 

V = P + (E + P)(p* + p* 2 P 2 ) = P + AP; (4) 

the density of mass-energy is 

£ = E + (E + P)(dE/dP)( P * +p* 2 P 2 ) = E + AE. (5) 

Here, pg and p 2 are dimensionless functions of R that are proportional to ft 2 , and describe the pressure perturbation, 
P is the pressure and E is the energy density of the non-rotating configuration. The stress-energy tensor for the fluid 
of the rotating object is 

r; = (f + %/-?^. (6) 

The rotational perturbations of the objects's structure are described by the functions u),ho,mo-,pQ,h 2 ,m 2 ,v 2 ,p 2 . 
These functions are calculated from Einstein's field equations (for details see fl7l fl8l|). 



B. The Exterior Solution 



The HT metric describing the exterior field of a slowly rotating slightly deformed object is given by 
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are the associated Legendre functions of the second kind. The constants M, J and Q are related to the total mass, 
angular momentum and mass quadrupole moment of the rotating object, respectively. This form of the metric corrects 
some misprints of the original paper by Hartle and Thorne [lH (see also J2f| and |2(| ) . 

The total mass of a rotating configuration is defined as M. = M' + SM, where M' is the mass of non-rotating 
configuration and SM is the change in mass of the rotating from the non-rotating configuration with the same central 
density. It should be stressed that in the terms involving J 2 and Q the total mass M. can be substituted by M' since 
SM is already a second order term in the angular velocity. 

In general, the HT metric represents an approximate vacuum solution, accurate to second order in the angular 
momentum J and to first order in the quadrupole parameter Q. In the case of ordinary stars, such as the Sun, 
considering the gravitational constant G and the speed of light c, the metric ([7]) can be further simplified due to the 
smallncss of the parameters: 
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For this special case one can calculate the corresponding approximate metric from ([7]) in the limit c 
computations are straightforward and lead to 



ds 2 



1 



2GM 



1 



c?R c 2 R 3 
2GM 2GQ 



2GQ , ^ 2G 2 MQ , n ' 
* Pi (cos 0) + 4p4 P 2 (cos 6) 



c 2 R c 2 R 3 



P 2 (cos6) 



dR 2 



c 4 i? 4 
_ 2GQ 
c 2 R 3 



c 2 dt 2 



AG J 
~^R, 



sin Qdtdtfi 



P 2 (cos9) 



(9) 

oo. The 
(10) 



R 2 (dQ 2 +sin 2 Od(f> 2 ) 



This metric describes the gravitational field for a wide range of compact objects, and only in the case of very dense 
(GM oj c 2 lZ) or very rapidly rotating (GJ ~ c 3 lZ 2 ) objects large discrepancies will appear. 



III. THE FOCK'S APPROACH 

Fock proposed in [20l | a method to analyze Einstein's equations in the presence of matter and to derive approximate 
interior and exterior solutions. This approach takes into account the internal properties of the gravitational source, 
and reduces the problem of finding interior approximate solutions to the computation of some integrals that depend 
explicitly on the physical characteristics of the object. In this section we present the main results of this approach 
and present a particular interior approximate solution and investigate the possibility of matching it with an exterior 
counterpart. 



The interior solution 



Fock's first-order approximation metric was recently re-derived and investigated by Abdildin [2 71 ] in a simple manner. 
Initially, this metric was written in its original form in a harmonic coordinate system [28l . [29| as follows 
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where U is the Newtonian gravitational potential, p is the mass density of the body, v is the speed of the particles 
inside the body, II is the comprcssional energy per unit mass, Pik is the stress tensor, U is the gravitational vector 
potential. The corresponding energy-momentum tensor is given as 
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Newton's potential satisfies the equation AU = —AirGp. The solution of this equation that satisfies the asymptoti- 
cally flatness condition at infinity can be written in the form of a volume integral: 
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Furthermore, the vector potential must satisfy the equation AUi = —AnGpvi whose general asymptotically fiat solution 
can be represented as 

Ui = G ( -j- l — dxidx' 2 dx' 3 . (15) 



Additional details about this metric can be found in [30| and [3l| . 

It is worth noticing that Chandrasekhar, using the Fock method, obtained in [32j a solution similar to (fTTj) that 
later on was used by Hartle and Sharp in |23|. However, it is not difficult to show that Chandrasekhar's solution is 
equivalent to (TlTl) . Indeed, the identification of the density 
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at the level of the energy-momentum tensor allows one to calculate the corresponding metric functions that show the 



equivalence of the metrics. Moreover, it has been shown in [32j that the solution for the non-rotating case can be 
matched with the well-known Schwarzschild solution, appropriately specialized to the case of spherical symmetry and 
hydrostatic equilibrium in the post Newtonian approximation. 

B. The exterior Fock metric for slow rotation and spherically symmetric distribution of mass 

In order to completely determine the metric, it is necessary to calculate the above integrals. Clearly, the result will 
depend on the internal structure of the body which is determined by the density p' and velocity v- distributions. Once 
these functions are give n, the calculation of the integrals can be performed in accordance with the detailed formalism 
developed by Fock [20| and then extended and continued by Abdildin [2ll, |22| and Brumberg [33|. Consider, for 
instance, the case of a rotating sphere with total mass M . Then, the corresponding exterior metric in spherical-like 
(non harmonic) coordinates can be written as plj 
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where So is the angular momentum of the body [53j . Here we added the constant k and verified that in fact the above 
metric is an approximate solution for any arbitrary real value of n. This simple observation allows us to interpret 
Fock's procedure as a method to find out how the internal structure of the object influences the values of the external 
parameters. For instance, the total mass in the above metric is M but it can decomposed as 

M = m+4, (18) 
cr 

where m is the rest mass of the body, and £ is an arbitrary real constant which, as the constant k, depends on the 
internal properties of the body. In particular, the cases of a liquid and a solid sphere have been analyzed in detail 
with the result 

l|T+|e, for a liquid sphere, _ I |, for a liquid sphere, 

1 4T + |e, for a solid sphere, 1 M' f° r a solid sphere. 

where T is the rotational kinetic energy of the body and e is the energy of the mutual gravitational attraction of the 
particles inside the body. 

C. The Kerr metric 

To describe the gravitational field of the rotating sphere outside the source, it seems physically reasonable to assume 
that the exterior vacuum metric be asymptotically flat. In this case, the first obvious candidate is the Kerr solution 
in the corresponding limit. The Kerr metric [1[ in Boyer-Lindquist coordinates [33l l3~ij can be written as 

ds 2 =U- 2 2 f 2q ) c 2 dt 2 - g l + f COs2l t d g 2 - (g 2 + a 2 cos 2 d) dd 2 

(20) 

2.2, 2/ipa 2 sin 2 \ . 2 2 A/jga sin 2 ■d 
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g z + a z cos z v J g z + a z cos z v 
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Expanding this metric to the order one obtains 
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Furthermore, if we introduce new coordinates g = g(r, 9), d = ?9(r, 9) by means of the equations 
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which coincides with the metric (fTT]) with k = 1. Consequently, the extended Fock metric (flTj) can be interpreted 
as describing the exterior field of a rotating body to second order in the angular velocity. The advantage of using 
Fock's method to derive this approximate solution is that it allows the determination of the arbitrary constant n. In 
fact, whereas k = kl = 4/7 for a liquid sphere and k = k$ — 15/28 for a solid sphere, the value for the Kerr metric 
k = kk = 1 does not seem to correspond to a concrete internal model. On the other hand, all the attempts to find a 
physically meaningful interior Kerr solution have been unsuccessful. Perhaps the relationship with Fock's formalism 
we have established here could shed some light into the structure of the interior counterpart of the Kerr metric. 
Furthermore, the coordinate transformation 11811 
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transforms the Kerr solution (|2T)|) . expanded to second order in the angular momentum, (here one should set 
G = c = 1) into the HT solution ([7]) with J = —A4a, AI = M. and a particular quadrupole parameter Q = J 2 j M.. 

In this way, we have shown that the extended Fock metric coincides for k = 1 with the approximate Kerr solu- 
tion which, in turn, is equivalent to the exterior HT solution with a particular value of the quadrupole parameter. 
The fact that in the Kerr solution the quadrupole moment is completely specified by the angular momentum is an 
indication that it can be applied only to describe the gravitational field of a particular class of compact objects. A 
physically meaningful generalization of the Kerr solution should include a set or arbitrary multipole moments which 
are not completely determined by the angular momentum. In the next section we present a particular exact solution 
characterized by an arbitrary quadrupole moment. 



IV. A SOLUTION WITH QUADRUPOLE MOMENT 



In this section we will consider the case of deformed objects as, for example, a rotating ellipsoid. It is obvious that 
if the form of the body slightly deviates from spherical symmetry, it acquires multipole moments, in particular, a 
quadrupole moment; the moments of higher order are negligible, especially, for a slowly rotating ellipsoid. We will 
generalize Fock's metric so that the quadrupole moment appears explicitly from the integration of pip and in the 
Newtonian potential. It should be mentioned that finding external and internal Newtonian potentials for a rotating 
ellipsoid is one of the classic problems of both physics and mathematical physics. Some examples for a homogeneous 
ellipsoid are considered in [35[ , but the most comprehensive details on this matter are given in (3(| and more recently 
in [371 ]. As for the exterior counterpart, there are several exact solutions [8HT^| with quadrupole moment and rotation 
parameter that could be used as possible candidates to be matched with the interior approximate solution. In this 
work, we limit ourselves to the study of a particular solution first proposed in [H] and then generalized in [l^, . 
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A. The exterior Fock solution 



Let us consider the first-order approximation metric pip . It is convenient to use the notation x\ = x, x 2 = y, and 
x 3 = z. In general, the fact that the mass density p — p(x, y, z) is a function of the coordinates does not allow one to 
find explicit expression for the internal Newtonian potential. It is possible only by numerical integration. However, 
for the case of uniform density there is in the literature a reasonable number of exact solutions for rotating ellipsoids. 
Since we consider slow rotation and the weak field approximation, we can use the expansion for the Newtonian 
potential (35|, [!(| 



U{r, 6)=G 



P' ,,,,,, Gm GD . 

■1—^rdx'dy'dz' = h — T P 2 cos( 

— r '\ r 2r J 



\r — r 



(27) 



where m is the rest mass of the ellipsoid, D is the Newtonian quadrupole moment, 6 is the angle between r = 
\J x 2 + y 2 + z 2 and z — axis. The first term in the expression above is the potential of a sphere and the second one is 
responsible for the deviation from spherical symmetry. If one takes the z axis as a rotating axis then the quadrupole 
moment is defined by 



D = J p{2z 2 - x 2 - y 2 )dxdydz. 



For the rotating ellipsoid with uniform density the quadrupole moment is 

D = ^(rl-rl) , 



(28) 
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where r p and r e are the polar and equatorial radii of the ellipsoid, respectively. The mass of the ellipsoid is defined 
as the integral m = J pdxdydz that in the case of an ellipsoid with uniform density yields 



m = -npr e r p 



(30) 



Note that the integration is carried out in the ranges of < x 1 y < r e and < z < r p . Using the same procedure one 
may write the integral in Fock's metric as follows 
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The quantity T>/c 2 is the relativistic correction to the Newtonian quadrupole moment D, i. e., the quadrupole moment 
due to rotation. 

To evaluate the integrals we use the relation for a compressible elastic medium [20j 



pU - pU + p = pil , 

where f2 is the potential of the centrifugal forces determined by 
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or v 2 = 2f2, which is a consequence of the equation Vi = UjiXj for the velocity within the body. In these equations 
ojik is the three-dimensional antisymmetric tensor of the angular velocity of the body. We may alternatively write its 
components as 



W23=0> a . = Wi, u 31 = LO y = uj 2 , Wia = uj z = w 3 . 
In our case 6j = {0, 0, tu z }; hence the potential of the centrifugal forces is 
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If we assume that the stress tensor within the body reduces to an isotropic pressure p, i. e., pik — ~pSik (a liquid 
always satisfies this condition), and take Eq. (j3~4"|) and v 2 = 2f2 into account, Eqs. (|32p and ([55)1 reduced to the simple 
form 

C = 2 / [2p£l + p] dxdydz , (38) 



V = 2 J [2pil + p] {2z 2 - x 2 - y 2 )dxdydz . (39) 

Furthermore, to evaluate these integrals we consider the following two cases that determine the inner structure of the 
body: 

1) A liquid body with following the equation of internal motion [20| 

f du \ dp fAn . 

''{07, 7T7,- (40) 

2) An absolute solid body with the following equation of internal motion [33| 

"T- = IT (41) 

OXi OXi 

In other words, these equations are simply the equations of hydrostatic equilibrium which are adopted by Fock and 
Brumberg to interpret the internal structure of an object. 

We limit ourselves to consider those cases in which the body rotates as a whole, in the manner of a rigid body. 
Then, for both liquid and solid bodies, the rotational kinetic energy takes the form 

T = J pfldxdydz = , (42) 

where I zz is the moment of inertia of the ellipsoid, which for a uniform density distribution, is equal to 

Izz = I p(x 2 + y 2 ) dxdydz = ^ mr l ■ (43) 



The pressure can be expressed as 



2T) , for a liquid body, 
for a solid body, 



where 



e = — J pU dxdydz , (45) 

represents the negative of the energy of mutual attraction of the constituent particles of the body. For a uniform 
density it has form 

3Gm 2 r p 

arccos — . (46) 
r e 



The second moments 

T ik = J pQxiX k dxdydz , (47) 
can be computed by using the above expressions. Then, for the second moments of the pressure we obtain (see [27| ) 

nx xh dxdydz - { ~ ^ Tlfc ' f ° r a Uquid b ° dy ' (48) 



2 



Tik, for a solid body, 



where (more details can be found in |27j ) 

2 

no. = — 5 



pxiXkXj — — dxdydz . (49) 
dxj 
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After calculating all the integrals we have 




for a liquid body, 
for a solid body, 



(50) 



V 



28 usSp 
5 ~7J7~ 



[/P(' 



y 2 )(z 2 
y 2 )(z 2 



x 2 )dxdydz\ — | J p(z 2 — x 2 )xj-§^-dxdydz 
x 2 )dxdydz] , 



for a liquid body, 
for a solid body, 



where So is the angular momentum of the body, which is found from 

Sq = I zz LO , 

and the numerical factors are 

for a liquid body, 
for a solid body. 

Unlike the Newtonian scalar potential, the vector potential can be easily calculated from 

G r. 




u 



2r 3 



So x f 



Introducing the effective (total) mass as 



M = m 

for the Fock metric we obtain the following expression 



(51) 



(52) 



(53) 



(54) 



(55) 



ds 2 



(GM GD 



P 2 (cos( 



\ r 2r 3 

2GM GD 
1+ — ^- + 



2 f GM GD 



P 2 (cos< 



[dr 



2r 3 
r 2 (d6 2 



P 2 (cos( 



P 2 (cos ( 



dt 2 



4GS 



(56) 



sin 8dcj)dt. 



This metric is written in harmonic coordinates. In order to write it in the Schwarzschild like (standard) spherical 
coordinates one should use the coordinate transformation 

GM 



which transforms the metric (|56p into 
, 2GM 



ds 2 



R 

2GM 



D + 
GD 



R 



%P2 (cos 6) 



e ->• e. 



c 2 R c 2 R 3 



P 2 (cos6) 



dR 2 



G DM 
c 2 R A 
GD 



c 2 R 3 



P 2 (cos0) 
P 2 (cos6) 



dt 2 



4:GS . 2AJ1M 

— 5 — sin <dd(pdt 
c z R 



(57) 



(58) 



P 2 (d6 2 +sin 2 9(# 2 ) , 



where we have neglected quadratic terms in the quadrupolc parameter D. In the limiting case with vanishing rotation 
Sq = and vanishing quadrupole moment D = T> = 0, this metric represents the approximate Schwarzschild solution. 

An examination of the metric (|56p shows that the rough approximation with r e m r p w r sp h ere and 5*0 ^ leads to 
the approximate Fock metric considered in Sec. |nl]with the total mass M , for a slowly rotating spherically symmetric 
body with 



D = 0, D 



2/^ (Sq 

M 



(59) 



It should be noted that an analogous result was obtained by Laarakkers and Poisson [4l[. They numerically 
computed the scalar quadrupole moment Q of rotating neutron stars for several equations of state (EoS). They found 
that for fixed gravitational mass M, the quadrupole moment is given as a simple quadratic fit 



Q 



-K 



J 



Mc 2 



(60) 



where J is the angular momentum of the star and K is a dimensionless quantity that depends on the EoS. The above 
quadratic fit reproduces Q with remarkable accuracy. The quantity K varies between K k. 2 for very soft EoS's and 
K « 8 for very stiff EoS's, for M = lAMsun as m neutron stars. This is considerably different from a Kerr black 
hole, for which K = 1 (see [Hj]). In our case, we have similar, but not the same results, since the Fock solution is not 
valid in the limit of strong gravitational fields (like in neutron stars) and fast rotation. The values for the constant 
k are obtained from qualitative analyses in the limit of weak field and slow rotation. In order to find exact values 
for k one should specify the EoS's and perform numerical integrations. This task, however, is out of the scope of the 
present work. 
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B. The exterior Quevedo-Mashhoon solution 

In this section, we study the general metric describing the gravitational field of a rotating deformed mass found by 
Quevedo and Mashhoon [l3l.l38l.l39l.l42j. which is a stationary axisymmetric solution of the vacuum Einstein's equations 
belonging to the class of Weyl-Lewis-Papapetrou [43M45| . For the sake of simplicity we consider here a particular 
solution involving only four parameters: the mass parameter M, the angular momentum parameter a, the quadrupolc 
parameter q, and the additional Zipoy-Voorhees constant 5. For brevity, in this section we use geometric units 

with G = c = 1. The corresponding line element in spheroidal coordinates (t, r, 9, <f) with r > a + Mo, < 9 < ^ is 
given by [48| 

ds 2 = f(dt ud<pf I je 2 ^ U 2 + r2 _ 2 ^ r + fl2 ) [(Mo - rf - (M 2 - a 2 ) cos 2 9} 

+ (r 2 ~2M r + a 2 ) sin 2 6>d0 2 1 , (61) 
where /, uj and 7 are functions of r and 9 only, and a is a constant. They have the form [x = (r — Mo)/a, y = cos 9] 

f = R-2 qS P 2 Q^ (62) 

u = -2a-2a^e 2 ^ 2 , (63) 

where 

R = a+a- +b + b_, L = a 2 + +b 2 + , (65) 

OH = (x + l) 5 " 1 [z(l - y 2 )(X + r))a+ + y{x 2 - 1)(1 - A ?/ )&+] , (66) 

7 = J(l + 9) 2 In + 2 «( 1 - P 2)Qi + 9 2 (1 - P 2 )[(l + P2){Q\ - Ql) (67) 

L x y 

+ ^(x 2 - 1)(2Q 2 - 3xQiQ 2 + 3QoQ 2 - Q' 2 )\. (68) 

Here Pj(j/) and Qi(x) are Legendre polynomials of the first and second kind, respectively. Furthermore 

a± = (x±l) B - 1 [x(l-\ri)±(l + \ri)], (69) 
b± = (x±l) 5 - 1 [y(X + V ) T (X-v)}, (70) 



with 



A = a(x 2 - l) 1 - s (x + y) 2S ~ 2 e 2qSS + 1 (71) 
n = a(x 2 -l) 1 - 5 (x-y) 2S - 2 e^ 55 -, (72) 

S± = \ln { -^yl + ^(l-y 2 T X y) + lKl-y 2 )Ty(x 2 -l)]ln^. (73) 



Morover, a and a are constants defined as 

a-M 



a = 



a = VA/ 2 - a 2 . (74) 



The physical meaning of the parameters entering this metric can be investigated in an invariant manner by calculating 
the Geroch-Hansen |49|, moments: 

M 2k+1 = J 2k = 0, k = 0,1,2,..., (75) 
M = M + a{5- 1), (76) 

M 2 = -Ma 2 + —ocr 3 - — cr((J - 1) [45M 2 + 15Ma(6 - 1) - (30 + 2q + 105 - 55 2 )a 2 ] , (77) 
15 15 

Ji = Ma + 2aa{5- 1), (78) 

J 3 = -Ma 3 + — agcr 3 - — acr(<S - 1) [60M 2 + 45Mcr(<5 - 1) - 2cr 2 (15 + 2q + 105 - 5c5 2 )l . (79) 
15 15 

The vanishing of the odd gravitoelectric (M„) and even gravitomagnetic ( J n ) multipole moments is a consequence of 
the reflection symmetry with respect to the equatorial plane 9 = tt/2. Note that in the limiting case 5 = 1, Mq = M 
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is the total mass of the body, a represents the specific angular momentum, and q is related to the deviation from 
spherical symmetry. All higher multipole moments can be shown to depend only on the parameters M, a, and q. In 
general, we see that the Zipoy-Voorhees parameter is related to the quadrupole moment of the source. In fact, even 
in the limiting static case with a = and q = 0, the only non-vanishing parameters are M = a and S so that all 
gravitomagnetic multipoles vanish and one obtains Mq = MS and Mi = — ^M 3 5(S 2 — 1) — the quadrupole moment 
that indicates a deviation from spherical symmetry. Some geo metrical properties of (|6ip versus particle motion and 
tidal indicators in this spacetime were explored in [25j | and 51 1, respectively. 

Consider the limiting cases of the QM solution. For vanishing quadrupole parameter, q = 0, 5 = 1, and vanishing 
angular momentum a = 0, a = 0, and a = M, one recovers the Schwarzschild solution with the following metric 
functions: 



2M 1 r(r- 2M) 

r 2 (M — r) 



)2 _ M 2 cos 2 , 



For vanishing quadrupole parameter and 5=1, one recovers the Kerr solution (|20p with i? 
functions 



(80) 

and g — ► r and 



2Mr 



2aMr$m 2 



1 r(r - 2M) + a 2 cos 2 6 

7 ~ 2 n (M-r) 2 -(A/ 2 -a 2 )cos 2 6l 



(81) 



The above limiting cases show that this solution describes the exact exterior field a rotating deformed object. To 
compute the case of a slowly rotating and slightly deformed body we choose the Zipoy-Voorhees parameter as 5 = l+sq, 
where s is a real constant. Then, expanding the metric (|ST|) to first order in the quadrupole parameter q and to second 
order in the rotation parameter a, we obtain 



2M 
r 



2a 2 M cos 2 I 



1 \2M 



- 1 



+ q{l + s) f 1 ) 

(3cos 2 0- 

r 



V2A7 



M 

1) (3cos 2 6»- 1) sin 2 

r 



In 1 - 



2M 



(82) 



2aA/rsin 2 6» 
r - 2Af ' 



(83) 



_ 1 r(r- 2M) 

1 ~ 2 n (r - M) 2 - M 2 cos 2 6 



a 



M 2 cos 2 sin 2 



+ + s)ln 



r(r - 2M)((r - A/) 2 - M 2 cos 2 0)_ 
r(r- 2 M) 



3q 



1 / r 
1+ 2(m 



The further simplification s 



r = R + Mq + -Mq sin 2 6 



2A4 5 



a 

2^ 



(r - 7\/) 2 - A/ 2 cos 2 
1, and the coordinate transformation fl8l. |25|. IB" 

R 2 



2M\ 



sin 2 0. 



(84) 



(85) 



= 9 - sin 6 cos 6 <; -q 



transforms the approximate QM solution 

ds 2 = 



into 

'2 



/ 2M \ 
2R\ + R ) 



2GM(l-q) 2G (J 2 4 
c 2 R c 2 R 3 \M 5 q 



AGMc 
c 2 R 



■ sin Odtdtfi 



1 + ^) ft(cose) 



3 2 rf^ 



1 - 



2G (3 



— --qM A P 2 (cos6) 



c 2 i? 3 IM 5 



i? 2 (d9 2 +sin 2 ed0 2 ) . 



(86) 



(87) 



Here we introduced again all the necessary constants G and c in order to compare our results with previous metrics. 
Finally, if we redefine the parameters M, a, and q as 



M = M(l-q), 3 = -Ma, Q = ^--^M 3 q, 



(88) 
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the approximate metric (|87|) coincides with the exterior HT metric (|10[) and, consequently, can be matched with the 
interior HT metric discussed in ITU 

The above metric is equivalent to the exterior extended Fock metric discussed in the previous subsection. To see this 
one has to consider the exterior solution ([55} which is written in the same coordinates as the exterior solutions (|57|) 
and (fTUf . It is convenient to show first the equivalence with the exterior HT solution (fTOj) that yields the conditions 



M = M, 



J — On 



1 



Q=-n\D+- S 



V 



The equivalence with the approximate QM solution (|87p follows then from the comparison of Eqs. 
obtain 



5 c 4 1 
8 G 2 AP 



D 



1 



V 



2J 2 
M 



and for vanishing D 



5 c 2 J 2 



(89) 

and (J5J]). We 

(90) 

(91) 



This result is in accordance with the limiting case of the Kerr metric for which we obtained that n = 1 and hence 
g = 0. 

Thus we come to the conclusion that in the limit of a slowly rotating and slightly deformed body the QM approximate 
solution is equivalent to the exterior Fock solution. 



V. CONCLUSIONS 



In this work, we studied the gravitational field of slowly rotating, slightly deformed astrophysical compact objects. 
We presented the main exact and approximate solutions of Einstein's equations that can be used to describe the 
interior and the exterior gravitational field. In particular, we presented the method proposed by Hartle and Thornc 
to find interior and exterior approximate solutions, and the method proposed by Fock to derive approximate interior 
and exterior solutions. We derived an extension of the approximate exterior Fock metric that takes into account up 
to the first order the contribution of a quadrupolc parameter that describes the deviation of the body from spherical 
symmetry. A particular parameter that enters the extended Fock metric turns out to have very specific values in the 
case of a liquid sphere and a solid sphere. In the case of the approximate Kerr metric, this parameter does not seem 
to correspond to any known interior model analyzed in the framework of Fock's formalism. 

We found that a particular QM solution, which in general possesses an infinite set of gravitational and electromag- 
netic multipolc moments, contains the exact Kerr metric and the approximate HT metric as special cases. Moreover, 
since the HT solution is endowed with its interior counterpart, we conclude that the approximate QM solution (to the 
second order in the angular momentum and to the first order in the quadrupole parameter) can be matched with the 
interior HT solution, indicating that it can be used to correctly describe the gravitational field of astrophysical com- 
pact objects. Moreover, we showed that the explicit form of the exterior Fock metric is equivalent to the approximate 
exterior QM solution. 

To avoid the technical problems that are usually found in the process of matching solutions, we use the same set 
of coordinates inside and outside the body. In the cases presented here, this can be done in a relative easy way 
only because all the coordinate transformations are not calculated exactly, but with the same approximation as the 
metric functions. This approach allows us to reduce the matching problem to the comparison of the metrics on the 
matching surfaces in such a way that only algebraic conditions appear. Using this method, we could show that the 
approximate Kerr metric cannot be matched with an interior Fock solution. However, if we take into account an 
additional quadrupole parameter, the matching of the extended Fock metric can be carried out by using as exterior 
counterpart a particular approximate QM solution that contains the Kerr metric as special case. We conclude that 
the quadrupole parameter offers an additional degree of freedom that allows the matching. It would be interesting to 
see if this is also true in the case of exact solutions. This could shed some light into the problem of finding a realistic 
gravitational source for the Kerr metric, a long-standing problem of classical general relativity. 
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